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We have performed an ab initio total energy investigation, within the density functional theory 
(DFT), of antisite defects in InP nanowires (InPiVlFs) grown along the [111] direction. Our total 
energy results indicate that, (i) P antisites (Pin) are the most likely antisite defect compared with 
In antisites (Inp), and (ii) the formation energies of P and In antisites do not depend on the NW 
diameter. In particular, thin InPiVWs, diameters of ~13 A, the Pi n antisite exhibits a trigonal 
symmetry, lying at 0.15 A from the Td site, followed by a metastable configuration with Pi n in an 
interstitial position (1.15 A from the Td site). We find a Pi n -P dissociation energy of 0.33 eV, and 
there is no EL2-like center for such a thin InPiVW. However, EL2-like defects occur by increasing 
the NW diameter. For diameters of ~18 A, the Pi n -P dissociation energy increases to 0.53 eV, 
which is 0.34 eV lower compared with the Pi n -P dissociation energy for InP bulk phase, 0.87 eV. 
We mapped the atomic displacements and calculated the relaxation energy, Franck-Condon shift, 
upon single excitation of Pi n induced states in InPiVIt 7 . The formation (or not) of EL2-like defects, 
Pin-P dissociation energy barrier, and the Franck-Condon (FC) energy shift, can be tuned by the 
NW diameter. 

PACS numbers: 71.15.Mb, 71.15.Nc, 71.20.-b 
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I. INTRODUCTION 

Antisite is one of the most studied native defect in III— 
V compounds, for instance, the EL2 center in GaAs. In 
a seminal work, Kamihska et al. [IJ established that the 
EL2 center exhibits a tetrahedral symmetry, formed by 
an isolated As antisite defect (Asc a )- Afterward, several 
experimental 0, Q as well as theoretical studies 0, H, 
@, H @| have been done aiming to clarify the structural 
and electronic properties of EL2 defect in GaAs. After 
more than ten year of investigations, currently there is a 
general agreement that the stable EL2 center is ruled by 
an isolated Asc a with Td symmetry, while the metastable 
EL2 M structure is attributed to the Asca atom displaced 
along the C$ v axis. 

Its is well known that low temperature growth of other 
III-V semiconductors, not only GaAs, allows the prepa- 
ration of highly nonstoichiometric compounds. Indeed, 
anion antisite defects in InP, Pi n , has been identified for 
the first time unambiguously through the electron param- 
agnetic resonance technique by Kennedy and Wilsey Q . 
Dreszer et al. using Hall, high-pressure far-infrared 
absorption and optically detected magnetic resonance 
measurements (at low temperature), verified that InP has 
two dominant donor levels associated with the phospho- 
rus antisite defect. Further investigations, based upon ab 
initio total energy calculations, proposed the formation 
of Pi n antisites clusters in InP bulk [ll| . 

The most interesting property on this class of defects is 
the fact that they also exhibit a clear metastable behav- 
ior [1, [12, [H[ • A neutral anion-antisite defect in III-V 
compounds has a stable fourfold and a metastable three- 
fold interstitial configuration, where the anion-antisite 



is displaced along the C$ v axis [l3|, similarly to EL2 
and EL2 M defects observed in GaAs. This center can be 
photo-excited into a metastable state, and from which it 
returns to the ground state by thermal activation. Thus, 
it is expected a persistent photoconductivity for Pi n at 
low temperature. 

Nowadays the electronic properties of several materi- 
als can be tailored throughout the manipulation/control 
of their size in the atomic scale. Within this new 
class of (nano)structured materials, the semiconductors 
nanowires (NWs) are attracting great deal of interest 
for future applications in several types of nanodevices. 
In particular, InP-nanowires (InPAVFs) have been con- 
sidered as a potential structure for fabrication of sensors, 
light emitting diodes, and field effect transistors [13, dl- 
Usually the vapor-liquid-solid mechanism, with a gold 
particle seed, has been utilized for the growth of these 
nanostructures (l6{ . These materials are quasi-one- 
dimensional with electrons confined perpendicularly to 
the NW growth direction. They exhibit interesting 
electronic and optical properties due to quantum con- 
finement effects, viz.: the size dependence of InPNW 
bandgap fl7l [Ta |. Recent ab initio calculations, per- 
formed by Li et al. , suggested the formation of stable 
DX center in small GaAs quantum-dots, dot diameter 
smaller than ~15 nm. 

It is quite likely the formation of native defects during 
the NW growth process. Therefore, the structural and 
electronic properties of those defects are important issues 
to be addressed, in order to improve our understanding of 
native defects in quasi-lD semiconducting NW systems. 

In this paper we carried out an ab initio total energy in- 
vestigation of antisite defects in InPNWs. We find that 
the formation of EL2-like defects, antisite dissociation 



energy, and the Franck-Condon (FC) energy are ruled 
by the NW diameter. On the other hand, the antisite 
formation energies are insensitive to the NW diameter, 
being Pi n antisites the most likely defect compared with 
Inp. At the equilibrium geometry, the Pi n atom in thin 
InPNW (diamenter of 13 A) exhibits an energetically 
stable trigonal symmetry, followed by a metastable con- 
figuration with Pi n in an interstitial position, 1.15 A from 
the Td site. In this case, the Pi n -P dissociation energy is 
0.33 eV. Increasing the NW diameter (13 — ► 18 A), we 
find a Pin-P dissociation energy of 0.53 eV, where Pi n ly- 
ing on the Td site represents the energetically most stable 
configuration. For the InP bulk phase, the Pi n -P disso- 
ciation energy is equal to 0.87 eV. Similarly to the bulk 
phase, Pin in InPNW induces the formation of localized 
states within the energy bandgap. However, for such a 
thin NW system, this defect does not exhibit an EL2-like 
behavior. On the other hand, increasing the NW diam- 
eter, EL2-like defects are expected to occur. Finally, 
within a constrained density functional approach [20j ]. 
we map the atomic displacements along thin InPNW 
upon single excitation of Pi n induced states, and calcu- 
late the respective relaxation energy, "Franck-Condon" 
(FC) shift. 
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FIG. 1: Structural models of thin InPAW, diamenter of 13 A. 
(a) Top view and (b) side view. The atomic geometry around 
the antisite defects are indicated in (c) Pi„, and (d) Inp. 



II. THEORETICAL APPROACH 

Our calculations were performed in the framework of 
the density functional theory (DFT) [2l[ , within the gen- 
eralized gradient approximation due to Perdew, Burke, 
and Ernzerhof 22). The electron-ion interaction was 
treated by using norm-conserving, ab initio, fully sepa- 
rable pseudopotentials [23[ . The Kohn-Sham wave func- 
tions were expanded in a combination of pseudoatomic 
numerical orbitals [24|. Double zeta basis set includ- 
ing polarization functions (DZP) was employed to de- 
scribe the valence electrons [25j. The self-consistent to- 
tal c harg e density was obtained by using the SIESTA 
code [2a]. The InPNW was modeled within the super- 
cell approach, where the InP bilaycrs were piled up along 
the [111] direction with periodicity length of 2\/3o and 
diameters of 13 and 18 A (a represents the optimized lat- 
tice constant along the [111] direction of InPNW). The 
NW surface dangling bonds were saturated with hydro- 
gen atoms. A mesh cutoff of 170 Ry was used for the 
reciprocal-space expansion of the total charge density, 
and the Brillouin zone was sampled by using one spe- 
cial k point. We have verified the convergence of our re- 
sults with respect to the number and choice of the special 
k points. All atoms of the nanowire were fully relaxed 
within a force convergence criterion of 20 me v/A. 

III. RESULTS AND COMMENTS 

Figured] presents the atomic structure of thin InPNW, 
diameter of 13 A, growth along the [111] direction, top 



view [Fig. QJa)] and side view [Fig. [ljb)] (the hydrogen 
atoms are not shown). Due to the ID quantum confine- 
ment, perpendicularly to the NW growth direction, the 
energy gap of InP NW increases compared with the bulk 
InP (1.0 eV) [H]. We find energy gaps of 2.8 and 2.2 eV 
for NW diameters of 13 and 18 A, respectively. It is 
important to take into account those energy gaps are un- 
derestimated, with respect to the experimental measure- 
ments, within the DFT approach. Figure [2ta) presents 
electronic band structure of thin InPAW for wave vec- 
tors parallel to the [111] direction (rL direction). At the 
r point, the highest occupied states exhibit an energy 
split of 0.31 and 0.20 eV (ai + e), for NW diameter of 
13 and 18 A, respectively. The valence band maximum 
of bulk InP is described by a three-fold degenerated t?, 
state. Such (£2 — > &i + e) energy splitting is due to the 
Td — > C^ v symmetry lowering of InP NW with respect to 
the bulk phase. 

The energetic stability of antisites in InPNWs can be 
examined by the calculation of formation energies (Oj). 
The formation energy of P and In antisites, Pi n [Fig.QJc)] 
and Inp [Fig.QJd)], respectively, can be written as, 

tti = E[InPNW l ) - E[lnPNW] - nmMin - np/ip. 

We have considered the formation of antisites in- 
ner InPNWs. E[lnP NWi] represents the total energy 
of InPNW with i = P\ n or Inp antisite defect, and 
E\hiPNW] is the total energy of a perfect InPNW. ni n 
(np) denotes the number of In (P) atoms in excess or in 
deficiency. The In and P chemical potentials, fj,i n and /ip, 
respectively, are constrained by following thermodynamic 
equilibrium condition, /ii n + /ip = /i^, P ' fc , where /i^ip is 
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FIG. 2: Electronic band structures of (a) perfect thin 
InPNW, and (b) defective thin InPNW with a P antisite 
(Pi„). NW diameter of 13 A. 



the chemical potential of bulk InP. Under In rich (P poor) 
condition we will have fii B — > nf^ lk , whereas under In 
poor (P rich) condition, /i In -> fJf™ lk - Af//(InP) 
For the heat of formation of bulk InP, Ai7/(InP), we 
have considered its experimental value, 0.92 eV. 

Figure [3] presents our calculated results of fii for thin 
InPNW (diameter of 13 A), as a function of the In chem- 
ical potential. It is clear that the formation of Pi n is 
dominant compared with Inp. This latter defect occurs 
only for In rich condition. At the In and P stoichiomet- 
ric condition (dashed line in Fig. [3]) we obtained f2p In = 
2.15 eV and £li np = 3.57 eV. Increasing the NW diameter 
(18 A), we find f2p In = 2.18 eV. Within the same calcu- 
lation procedure, we obtained similar formation energy 
results for the InP bulk phase, viz.: fip In = 2.14 eV and 
^in P =3.59 eV. Those results indicates that, (i) there is 
an energetic preference of Pi n antisites, compared with 
Inp , for both structural phases of InP. (ii) The formation 
energy of Pi n does not depend on the NW diameter. 

At the equilibrium geometry, Pi n defect in bulk InP 
keeps the Td symmetry with Pi n -P bond lengths of 
2.49 A, while Inp exhibits a weak Jahn-Teller distor- 
tion along the [001] direction. Those results are in accor- 
dance with previous ab initio studies of antisites in InP 
[ll],[28|, 29] . On the other hand, the equilibrium geometry 
of antisites in thin InPiVW is quite different. The Pi n -Pi 
bond, parallel to the NW growth direction, is stretched 
by 27% compared with the other three Pi n — P« bonds, 
i = 2 — 4 in Fig. [IJc). Similarly for the Inp antisites, 
Inp-Ini is stretched by 8.6% compared with the other 
three Inp-In bonds indicated in Fig. QJd). Figures QJa) 
and 0Jb) depict the total charge densities along the Pi n - 
Pi and Inp-Ini bonds, respectively. In particular, due to 
the large Pi n -Pi bond stretching, Pr n is weakly bonded 
to Pi, whereas the other Pi n -P2, -P3 and -P4 bonds [see 
Fig. Ufa)] exhibit a strong covalent character. Therefore, 
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FIG. 3: Formation energies of Pi n and Inp antisites in thin 
InPiVW, diameter of 13 A, as a function of In chemical po- 
tential. 



different from bulk InP, Pj n antisites in NW system ex- 
hibits a 631, symmetry, with the Pt„ atom displaced from 
Trf site by 0.15 A along the [111] axis. Such a Pi n dis- 
placement, along the C^ v axis, has not been observed by 
increasing the NW diameter to 18 A. In this case, the 
Pin atom accupies a Td site and the Pi n ~P bond lengths 
are equal to 2.94 A. 




FIG. 4: Total charge densities of (a) Pi n , and (b) Inp antisites 
in InPiVW. The partial charge densities of Pin induced (c) 
highest occupied vO and (d) lowest unoccupied cl states. 

Several anion antisite defects in III- V materials were 
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studied by Caldas et al. [HI]. In that work, based 
upon ab initio total energy calculations, the authors ob- 
served an energetically stable Td symmetry for Pj n and 
a metastable C% v configuration for the antisite atom dis- 
placed by ^1.3 A along the [111] direction. Indeed, in 
Fig. 03a) we present our total energy results as a func- 
tion of Pi n displacement, along the [111] direction, for 
the InP bulk phase. We find that the Td symmetry rep- 
resents the energetically most stable configuration, fol- 
lowed by an energy barrier of 0.87 eV at 0.7 A from the 
Td site, z = 0.7 A, breaking the Pi n -Pi bond. Finally, the 
metastable C^ v configuration occurs for z = 1.2 A, where 
the Pi n occupies an interstitial site. Figures Ob) and 
[5tc) present the energy barrier for Pi n in thin InPNW. 
We have examined two different processes: (i) the In and 
P atoms of the NW are not allowed to relax during the 
Pin displacement [Fig. GD[b)]. (ii) The In and P coordi- 
nates are fully relaxed for each Pi„ step along the [111] 
direction, i.e. an adiabatic process [Fig. G2c)]. It is no- 
ticeable that the energy barrier calculated in (i) is very 
similar to that obtained for Pi n in bulk InP [Fig. OJa)]. 
In (i) the energetically most stable configuration for Pi n 
exhibits a Td symmetry, and we find a dissociation energy 
of 0.84 eV for z = 0.86 A. Further Pj n displacement in- 
dicates a Cs v metastable geometry for z w 1.15 A, where 
the Pin atom occupies an interstitial site. Meanwhile, in 
(ii) the C^ v symmetry, with the Pi n atom lying at 0.15 A 
from the Td site, represents the energetically most stable 
configuration. The Pi n -Pi dissociation energy reduces to 
0.33 eV (at z = 0.85 A), and there is a metastable geom- 
etry for z w 1.15 A. The total energy difference between 
the stable (S) and the metastable (M) configurations, 
E{S) - E(M), in InPNW is equal to -0.084 eV, while 
in bulk InP we find E(S) - E(M) = -0.53 eV. Com- 
paring (i) and (ii) we verify that the atomic relaxation 
plays a fundamental rule not only to the energy barrier 
for the Pin atomic displacement along the [111] direction, 
but also for the equilibrium geometry of the Pi n struc- 
ture. Increasing the NW diameter, the energy barrier for 
Pin displacement will become similar to that obtained for 
bulk InP (an upper limit for very large NW diameter), 
see Fig.[5][d). Indeed, for NW diameter of 18 A, we find a 
Pin-Pi dissociation energy of 0.53 eV, and E(S)-E(M) 
equal to —0.14 eV. Those results indicate that the Pi n -Pi 
dissociation barrier, and the E(S) — E(M) total energy 
difference, can be tuned within the shaded region indi- 
cated in Fig. E3d). 

Focusing on the electronic structure, the formation of 
Pi n defect in thin InPNW gives rise to a deep double 
donor state, labeled vO in Fig.^b), lying at 1.3 eV above 
vl. Such a Pi n induced state is very localized within the 
fundamental band gap, being almost flat along the TL di- 
rection. For Pin in bulk InP, ab initio studies performed 
by Seitsonen et al. [28] indicates the formation of an oc- 
cupied state at 0.7 eV above the valence band maximum. 
The partial charge density contour plot of vO, Fig. BJc), 
shows an anti-bonding n* orbital concentrated along the 
Pin-Pi bond. The lowest unoccupied state, cl, is also 



concentrated on the Pi„ antisite, as depicted in Fig.QJd). 
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FIG. 5: Total Energy as a function of the Pi n antisite dis- 
placement along the [111] direction: (a) bulk InP, and (b) 
InPiVW where the In and P atoms are not allowed to relax 
during the Pi n displacement, (c) Whole In and P atomic po- 
sition along the NW are allowed to relax. In (c) the total 
energy barrier indicated by triangles (dashed line) was calcu- 
lated for an excited electronic configuration, single excitation, 
(d) Energy barriers for bulk InP and thin InPNW systems. 

EL2 and EL2 M centers in III-V are characterized by a 
stable Td and a metastable C$ v geometries for an isolated 
Vin interstitial atom. However, before the Vm antisite 
arrives to the metastable C^,v configuration, there is a 
strong electronic coupling between the highest occupied 
and the lowest unoccupied antisite induced states near 
to the local maximum for the VnrV dissociation energy, 
see Fig. 3 in Ref. [3 and Fig. 4 in Ref. 

Similarly, for thin InPNW , we examined the energy 
positions of single particle eigenvalues, e(v0) and e(cl), 
as a function of the Pi n displacement along the C% v axis. 
Our calculated results, depicted in Fig. ® reveal that 
there is no electronic coupling between vO and cl. Thus, 
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based upon the calculated energy barrier [Fig. [5[c)], and 
the evolution of the antisite induced states, vO and cl, 
we can state that there is no EL2-like center in thin 
InPNWs. On the other hand, increasing the NW diam- 
eter [Fig. [S^d)], EL2-like center is expected to occur in 
InP. 



> 



1.5 




X cl 




1.0 


r InPAW 






0.5 




relaxed 












0.0" 








-0.5 














■ ■ 











0.0 1.0 

P displacement (angs.) 



FIG. 6: Single particle eigenvalues (e) for Pi n : highest occu- 
pied vO [Fig. \M. C )] an d lowest unoccupied cl [Fig. |4jd)]. 

Within the DFT approach, we have examined the 
structural relaxations upon single excitation from the 
highest occupied state (vO) to the lowest unoccupied state 
(cl) of Pin in thin InPNW. We have used the calcula- 
tion procedure proposed by Artacho et al., where the 
single excitation was modeled by "promoting an electron 
from the highest occupied molecular orbital (HOMO) to 
the lowest unoccupied molecular orbital (LUMP)" [2fj|. 
i.e. a "constrained" DFT calculation. So that, promot- 
ing one electron from vO to cl, vO — > cl, we obtained 
the equilibrium geometry (full relaxations) as well as the 
self-consistent electronic charge density. The atomic re- 
laxations along thin InPNW, upon vO — > cl single exci- 
tation, are localized near to the Pi n position. Since we 
are removing one electron from the anti-bonding n* or- 
bital depicted in Fig. BJc), the Pi n -Pi bond shrinks by 
0.12 A (2.90 — > 2.78 A), and there is radial a contraction 
of -0.2 A at z = 6.8 A. 

The structural relaxations are proportional to the 
degradation of optical energy due to the atomic displace- 
ments associated with such a vO — > cl single excitation, 
FC shift. Comparing the radial and the longitudinal dis- 



placements, we verify that the most of the FC shift comes 
from the atomic relaxation parallel to the NW growth 
direction. The atomic displacements along the radial di- 
rection (Ar) are localized nearby Pi n , lying within an 
interval of —0.2 < Ar < A. While the atomic displace- 
ments parallel to the growth direction (Az) lie within a 
range of -0.2 < Az < 0.3 A, being less concentrated 
around the Pi n antisite position. We find a FC energy 
shift of 0.8 eV, which is quite large compared with the 
InP bulk phase. Within the same calculation approach, 
we find a FC energy shift of 0.05 eV for InP bulk. Previ- 
ous ab initio study indicates a FC energy shift of 0.1 eV 
for bulk InP. Thus, suggesting that the FC energy shift 
also can be tuned by controlling the NW diameter. 

IV. CONCLUSIONS 

In summary, we have performed an ab initio total en- 
ergy investigation of antisites in InPNW. We find that 
the Pi n antisite is the most likely to occur. For thin 
InP TV Ws the Pi„ atom exhibits a trigonal symmetry, fol- 
lowed by a metastable P interstitial configuration. In- 
creasing the NW diameter (18 A), Pi n occupying the Td 
site becomes the energetically most stable configuration. 
The calculated energy barrier, for Pi n displacement along 
the C% v axis, indicates that there is no EL2-like defect 
in thin InPNWs. However, increasing the NW diame- 
ter, 13 — > 18 A, we observe the formation of EL2-like 
defects in InP. Within a "constrained" DFT approach, 
we calculated the structural relaxations and the FC en- 
ergy shift upon single excitation of the electronic states 
induced by Pi„. We inferred that not only the formation 
(or not) of EL2-like defects, but also (i) the Pi n -P dis- 
sociation energy, (ii) the total energy difference between 
stable and metastable Pi n configurations, E(S) — E(M), 
and the (iii) FC energy shift can be tuned by controlling 
the InP AW diameter. 
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